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ABSTRACT 

We report the detection of GeV 7-ray emission from the supernova remnant Puppis A 
with the Fermi Gamma-Ray Space Telescope. Puppis A is among the faintest supernova 
remnants yet detected at GeV energies, with a luminosity of only 2.7xl0 34 (D/ 2.2 
kpc) 2 erg s -1 between 1 and 100 GeV. The 7-ray emission from the remnant is spatially 
extended, with a morphology matching that of the radio and X-ray emission, and is 
well described by a simple power law with an index of 2.1. We attempt to model the 
broadband spectral energy distribution, from radio to 7-rays, using standard nonthermal 
emission mechanisms. To constrain the relativistic electron population we use 7 years of 
WMAP data to extend the radio spectrum up to 93 GHz. Both leptonic and hadronic 
dominated models can reproduce the nonthermal spectral energy distribution, requiring 
a total content of cosmic ray (CR) electrons and protons accelerated in Puppis A of at 
least Wcr ~ (1 - 5)xl0 49 erg. 
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1. Introduction 

Supernovae have long been thought responsible for accelerating protons to relativistic energies 
in our Galaxy. Diffusive shock acceleration provides a mechanism to convert a significant fraction 
of the blast wave kinetic energy into relativistic particles. However it has proven difficult to mea- 
sure cosmic rays (CRs) from supernova remnants (SNRs) directly (Reynolds 2008, and references 
therein) . 

The Fermi Gamma-Ray Space Telescope (Fermi) has identified a number of SNRs as high- 
energy y-ray sources. Middle-aged SNRs interacting with molecular clouds are among the most 
luminous GeV sources in the Galaxy (Abdo et al. 2009a, 2010a, e,f; Castro & Slane 2010). The 
high luminosities are thought to result from interactions between cosmic rays and the large target 
mass of the molecular clouds, though it is not yet certain whether GeV emission arises from lep- 
tonic (Bremsstrahlung) or hadronic (7r° decay) processes. In contrast, the young SNR Cassiopeia 
A (hereafter Cas A), is detected with a luminosity two orders of magnitude lower than SNRs in- 
teracting with molecular clouds (Abdo et al. 2010c). A careful consideration of y-ray emission 
mechanisms in Cas A indicates that at least 6xl0 49 erg in hadronic CRs are required (Araya & Cui 
2010). Other young SNRs have also been detected by Fermi - Tycho (Giordano et al. 2012), RX 
J1713.7— 3946 (Abdo et al. 2011), and RX J0852. 0-4622 (Tanaka et al. 2011) - though the origin 
of y-ray emission remains uncertain. Clearly, SNRs contribute at least part of the Galactic CR 
population, though the total energetic input and timescale for acceleration require further study. 

Supernova remnant Puppis A (G260.4-3.4) is an important case to study, as it shows signs 
of recently encountering a higher ambient density in the vicinity of a nearby molecular cloud 
(Hwang et al. 2005). The remnant displays increasing X-ray surface brightness from west to east 
(Petre et al. 1982) corresponding to an increasing density of the ambient interstellar medium (ISM) 
at the eastern and northern shell (Dubner & Arnal 1988). The proper motions of fast optical 
knots gives a dynamical age of 3700±300 years, establishing that the SNR is in the Sedov-Taylor 
evolutionary phase (Winkler et al. 1988). X-ray spectra from the rim of Puppis A show low 
elemental abundances relative to solar, consistent with swept-up interstellar medium (ISM) and a 
lack of ejecta contamination (Katsuda et al. 2008). A progenitor mass of 15-25 M 0 is inferred from 
comparisons of relative abundances of metal-rich ejecta with nucleosynthesis models (Hwang et al. 
2008; Katsuda et al. 2010). A plausible SN IIL/b progenitor would have produced a clumpy, red 
supergiant wind only out to 7 pc (Chevalier 2005). At an estimated distance of 2 kpc (Reynoso 
et al. 1995), the diameter of Puppis A is 30 pc, so the remnant is no longer interacting with the 
circumstellar medium of the progenitor, but the surrounding ISM. 

The extent of the shock interaction is characterized by thermal dust emission in infrared (IR) . 
Arendt et al. (1991) noted enhanced dust emission corresponding with the bright X-ray shell. The 
dust emission is approximated by a 65 K black body, whereas thermal X-rays originate from a 
>10 6 K gas (Hwang et al. 2005). Spitzer Space Telescope imaging of the SNR at 24, 70 and 160 
/am reveals a detailed correlation between IR and X-ray emission on arcsecond scales (Arendt et al. 
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2010). This is interpreted as swept-up interstellar dust heated by collisions with the hot shocked 
plasma observed in X-rays. There are two notable regions where the SNR shock has engulfed small 
denser clouds: the Bright Eastern Knot and the Northern Knot. Besides those the shock has not 
yet become radiative for most of the SNR, consistent with the relatively young age and the low 
density of the surrounding medium. 

Here we report the detection of high energy 7 -ray emission from 200 MeV to 100 GeV, based 
on observations with the Large Area Telescope (LAT), the primary science instrument on Fermi. 
We also analyze the WMAP 7-year data from 23 to 93 GHz, in order to better constrain the 
relativistic electron population. The excellent characterization of the physical properties of the 
Puppis A makes it an ideal candidate to explore the production of nonthermal emission this SNR. 


2. Observations and Data Analysis 
2.1. Fermi- LAT observations and data reduction 

The LAT detects 7 -ray photons by conversion into electron-positron pairs in the energy range 
between 20 MeV to higher than 300 GeV, as described by Atwood et al. (2009). It contains a high- 
resolution converter/tracker (for direction measurement of the incident 7 -rays), a CsI(Tl) crystal 
calorimeter (for energy measurement), and an anti-coincidence detector to identify the background 
of charged particles. The LAT has a large effective area (^ 8000 cm 2 on-axis above 1 GeV), a wide 
field of view (~ 2.4 sr) and good angular resolution (^0.6° radius for 68 % containment at 1 GeV 
for events converting in the front section of the tracker). The on-orbit calibration is described in 
Abdo et al. (2009b). 

The following analysis was performed using 36 months of data collected from 2008 August 4 
to 2011 August 20 within a 15° x 15° region around the position of Puppis A. Only events with 
Earth zenith angles smaller than 100° were included to reduce contamination from the Earth limb. 
We used the P7V6 instrument response functions (IRFs), and selected the ‘Source’ events which 
correspond to the best compromise between the number of selected photons and the charged particle 
residual background for the study of point-like or slightly extended sources. 

Two different tools were used to perform the spatial and spectral analysis: gtlike and 
pointlike, gtlike is a binned maximum-likelihood method (Mattox et al. 1996) implemented 
in the Science Tools distributed by the Fermi Science Support Center (FSSC) 1 . pointlike is an 
alternate binned likelihood technique, optimized for characterizing the extension of a source (unlike 
gtlike), that was extensively tested against gtlike (Kerr 2011; Lande et al. 2012). These tools fit 
a source model to the data along with models for the residual charged particles and diffuse 7 -ray 


^^More information about the performance of the LAT can be found at the FSSC 
(http : //f ermi.gsfc.nasa.gov/ssc). 
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emission. In the following analysis, the Galactic diffuse emission is modeled by the standard LAT 
diffuse emission ring— hybrid model gaL2yearp7v6-v0.fits. The residual background and extragalac- 
tic radiation are described by a single isotropic component with a spectral shape described by the 
file iso_p7v6source.txt. These models are available from the FSSC. 

Sources within 15° around Puppis A listed in the Fermi - LAT Second Source Catalog (Nolan 
et al. 2012, hereafter 2FGL) are included in our spectral-spatial model of the region. We also 
include extended spatial templates describing the y-ray emission from RX J0852. 0-4622 (Tanaka 
et al. 2011) and the pulsar wind nebula (PWN) Vela-X (Grondin et al. 2012). Puppis A is associated 
with three 2FGL sources within the radio boundary: 2FGL J0821. 0-4254, J0823. 0-4246, J0823.4- 
4305. These sources are either all left free to vary or are replaced by other spatial models. The 
spectral parameters of sources closer than 5° to Puppis A are left free, while the parameters of all 
other sources are fixed at the 2FGL values. 

Located at a distance of 3° from the SNR Puppis A, the Vela pulsar is the brightest steady 
y-ray source in the sky, from which photons are observed up to 25 GeV (Abdo et al. 2010g). To 
avoid any bias on the analysis of the SNR due to this bright nearby pulsar, all studies involving 
events below 3 GeV were performed in the off-pulse window of the Vela pulsar. Gamma-ray photons 
were phase-folded using a timing solution produced from observations made with the Parkes 64 m 
radio telescope (Weltevrede et al. 2010). To build this timing solution, 163 times of arrival (TOA) 
were used covering the period from the beginning of the mission to 2011 August 20. We fit the 
radio TOAs to the pulsar rotation frequency and first ten derivatives. The fit further includes 15 
harmonically related sinusoids, using the FITWAVES option in the TEMP02 package (Hobbs et al. 
2006), to flatten the timing noise. The post-fit rms is 119.579 /is, or 0.13% of the pulsar phase. 
This timing solution will be made available through the FSSC 2 . 

The off-pulse phase window was chosen between 0.65 - 1.05, slightly larger than the one used 
for the analysis of the Vela-X PWN (Abdo et al. 2010b) allowed by the larger angular distance 
between Puppis A and Vela than between the pulsar and its associated PWN. This off-pulse choice 
offers good statistics without increasing the systematics due to the contamination from the Vela 
pulsar. 


2.2. Morphological analysis of Fermi- LAT data 

In the study of the morphology of an extended source, one of the major objectives is to have 
the best possible angular resolution. Therefore, we restrict our LAT dataset to events with energies 
above 800 MeV. This also reduces the relative level of the Galactic diffuse background owing to its 
softer spectrum. As discussed above, to ensure that we do not suffer contamination from the Vela 
pulsar, we only selected the off-pulse photons defined as having phases 0.65 - 1.05. 


2 Public page for pulsar timing models: http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
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We fit models to the data using the maximum likelihood framework. The test statistic (TS) 
of a source is defined as twice the logarithm of the ratio between the likelihood C\ obtained by 
fitting the source model plus background components (including other sources) to the data, and the 
likelihood obtained by fitting the background components only, i.e TS = 2 log(£i/£o)- Figure 1 
shows a map of the point-source detection significance, evaluated at each point in the TS map for 
the region around Puppis A using photons with energies > 800 MeV. This skymap contains the TS 
value for a point source of fixed photon index T = 2 at each map location, thus giving a measure 
of the statistical significance for the detection of a y-ray source with that spectrum in excess of 
the background. Emission coincident with Puppis A is detected, but also in two regions outside 
of its shell. We consider these regions of excess y-ray emission as likely to belong to background 
sources not recognized in 2FGL. Due to the longer integration time of our analysis (36 months vs. 
24 months in the catalog) and the overwhelming brightness of the Vela pulsar in the entire phase 
interval, the appearance of additional sources in our region of interest is expected. These sources 
are denoted with the identifiers BckgA and BckgB and will be described below. 

Puppis A was first established as an extended source above 1 GeV by Lande et al. (2012). 
Following their procedure, we determined the source extension using point like with a uniform 
disk hypothesis (compared to the point-source hypothesis) for energies above 800 MeV in the off- 
pulse of the Vela pulsar and above 5 GeV in the all-phase interval. The results are summarized 
in Table 1. To quantify the significance of the extension of the SNR Puppis A, we define TS ex t as 
twice the difference between the log-likelihood of an extended source model and the log-likelihood 
of a point-like source model. Above 800 MeV, TS e ^=46 (which converts into a significance of ^ 7a 
for the source extension), demonstrating that the source is significantly extended with respect to 
the LAT point-spread function (PSF). The fitted radius is 0.38° db 0.04°, in good agreement with 
the size of the SNR as seen in the radio and X-rays. The extension is still significant above 5 GeV 
with a slightly smaller size of 0.32° ± 0.03°. Going from lower to higher energy, the radius of the 
SNR seems to decrease with a fitted position centered on the Eastern side of the remnant but, due 
to the large uncertainty of the fit with the current statistics, this effect is only marginal. This trend 
is visible in Figure 2 which presents the LAT TS maps above 800 MeV and above 5 GeV when 
nearby sources and the extragalactic and Galactic components of the background are included in 
the background model. The emission above 5 GeV is brightest in the East, near the bright eastern 
knot, where the SNR shock is interacting with a dense cloud, and similar to the radio and X-ray 
morphologies. 

It should be noted that during the fit above 800 MeV, we simultaneously localized the sources 
BckgA and BckgB assuming the best disk spatial model for the SNR Puppis A, presented in 
Table 1. We followed the procedure described in Lande et al. (2012). BckgA was fit to position 
<a(J2000) = 125.77° , 5(J2000) = —42.17° with a 68% error radius of 0.06°, while BckgB was fit at 
<a(J2000) = 128.14°, 5(J2000) = —43.39° with a 68% error radius of 0.05°. The addition of sources 
BckgA and BckgB in the source model results in an improved maximum likelihood above 800 MeV 
with four additional degrees of freedom. The sources have corresponding significances of ^ 4.4 and 


- 6 - 


^ 5.5a (TS = 27 and 39) respectively, demonstrating that these two sources are significant and 
distinct from Puppis A. In addition, we will show in § 2.3 that, in a combined likelihood analysis 
of the spectrum of Puppis A, the emission from these two sources is much softer than the y-ray 
emission from the SNR. 

To be more quantitative, we have examined the correspondence of the y-ray emission from 
Puppis A with different source morphologies by using gtlike with assumed multi- frequency tem- 
plates above 800 MeV. For this exercise we compared the TS of the best-fit uniform disk model 
provided by point like (see Table 1) with TS derived when using morphological templates from 
ROSAT and the VLA (1.4 GHz radio images) assuming a power-law spectrum. The resulting TS 
values obtained from our maximum likelihood fitting are summarized in Table 2. The radio tem- 
plate, the X-ray template and the uniform disk all produce improvements from the 2FGL 3-point 
source model, while having fewer degrees of freedom. Finally, we divided the ROSAT and uniform 
disk templates into two half-disks along the North/South line in celestial coordinates, as indicated 
in Figure 1, to verify the evidence from point like of an energy-dependent shape. The two-half- 
disk model provides a marginal improvement (^ 2.4cr level), in the same trend described above. 
However, further data are needed to confirm this effect. 


2.3. Spectral analysis of Fermi- LAT data 

Using the different templates described above, we performed maximum likelihood fits and 
compared the best-fit parameters in the wide energy range 0.2-100 GeV. To avoid any bias due to 
the brightness of the Vela pulsar, the spectral fit was performed in the off-pulse phase interval. No 
evidence for cutoff or break is visible and, in all cases, the Fermi - LAT data are well described by 
a power-law such that the differential photon flux, F(E ), is given by: 


F(E) 


dN 

Je 



( 1 ) 


where Fq is the differential flux at energy Fq and T is the photon index. The uncertainty contours 
as a function of energy F, called the butterfly, are defined such that the differential flux satisfies: 


A F 2 
F 2 



+ ‘* 2 (I) Ar2 


2 

Fb 


cov(F 0 ,r) l n 
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where cov(Fq, T) is the covariance term, returned by the MINUIT minimization and error analysis 
function called by the Fermi-L AT likelihood analysis tool, gtlike, and A F, AFo and AY are the 
statistical uncertainties on F, Fq and T, respectively. 

As shown in Table 3, the obtained fluxes and spectral indexes are almost the same for all 
spatial templates used. These spectral fits take into account the y-ray emission of the two additional 
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Table 1. Centroid and extension fits to the LAT data for Puppis A using point like for events 
with energies between 800 MeV and 100 GeV and between 800 MeV and 5 GeV in the off-pulse 
window of the Vela pulsar, and above 5 GeV in the all-phase interval. 


Spatial Model 

Name 

Energy (GeV) 

R.A. (°) 

Dec. (°) 

Radius (°) 

TS e z£ 

Point Source 

PS 1 

0.8 - 100 

125.85 

-42.90 





Disk 

D 1 

0.8 - 100 

125.67 

-42.93 

0.38 

± 

0.04 

46 

Point Source 

PS 2 

0 . 8-5 

125.62 

-42.92 





Disk 

D 2 

0 . 8-5 

125.71 

-42.85 

0.47 

± 

0.08 

14 

Point Source 

PS 3 

5 - 100 

125.83 

-42.87 





Disk 

D 3 

5 - 100 

125.69 

-42.90 

0.32 

± 

0.03 

26 


Table 2. Test statistic obtained using gtlike for different spatial models compared with the null 
hypothesis of no 7 -ray emission associated with the SNR Puppis A. The photon energies are 0.8 - 
100 GeV; only off-pulse photons are selected. For each model, we give the name, the difference in 
log-likelihood value and the number of additional degrees of freedom (Ndf) 


Spatial Model 

TS 

Ndf 

Null hypothesis 

0 

0 

Point Source (PS1) 

120 

4 

Three Point Sources (2FGL) 

160 

12 

Radio Template 

166 

2 

X-ray Template 

170 

2 

Uniform Disk (Dl) 

172 

5 

Two Half-Disk (E, W) 

180 

7 

Two X-ray halves (E, W) 

178 

4 


background sources BckgA and BckgB. Their positions and spectral parameters (fitted using the 
uniform disk D 1 to describe Puppis A) are given in Table 4. 

We searched for spectral variations across the 7 -ray emission associated with Puppis A. We 
allowed an independent normalization and spectral index for the two X-ray halves. Interestingly, 
the spectrum of the Western half-disk is found to be steeper by an index of 0.3 (2a) above 200 
MeV, with respect to the Eastern half-disk. However, more statistics are needed to confirm spectral 
differences between the Eastern and Western regions. The spectral parameters for each of the two 
half-disks of the remnant are reported for reference in Table 3. Note that the TS of each individual 
half is computed by removing only that half, so the total likelihood improvement from the two 
X-ray halves model (Table 2 ) is more than the sum of the TS of the individual halves (Table 3). 

To ensure that the slightly softer spectrum in the Western half-disk is not due to contamina- 
tion from the compact object inside the SNR Puppis A, PSR J0821— 4300, we performed a timing 
analysis using the same Fermi - LAT dataset. The low TS of this side of the remnant combined 
with the need to keep only photons in the Vela pulsar off-pulse make a blind pulsation search unfea- 
sible. We used ephemerides from Gotthelf & Halpern (2009) with the assumption of a zero period 
derivative (Gotthelf et al. 2010) to fold the photon arrival times using the Fermi plugin distributed 
with the TEMP02 software. We applied the H - test periodicity test (de Jager &; Busching 2010) 
for several cuts in both energy and photon distance from the pulsar location and found no evidence 
of pulsation with a significance better than H- test=13.5 (2.8cr). This confirms that with a small 
spin-down luminosity and a relatively weak magnetic field (Gotthelf et al. 2010) PSR J0821— 4300 
is not a good candidate for 7 -ray emission. 

The Fermi - LAT spectral points for Puppis A were obtained by dividing the 200 MeV - 100 
GeV range into seven logarithmically-spaced energy bins and performing a maximum likelihood 
spectral analysis in each interval, assuming a power-law spectral shape with fixed photon index 
T = 2 and the X-ray spatial model for Puppis A. Spectral points of Puppis A above 3 GeV were 
determined using the all-phase interval which is allowed by the improved PSF at those energies 
as well as the lower energy flux of the Vela pulsar. These spectral energy points are in agreement 
with those derived in the off-pulse window. The result and the butterfly from the overall spectral 
fit, renormalized to the total phase interval, are presented in Figure 3 along with the butterflies 
corresponding to the fit of the two half-disks of the SNR. The errors on the spectral points represent 
the statistical and systematic uncertainties added in quadrature. Two main systematic errors have 
been taken into account: imperfect modeling of the Galactic diffuse emission and uncertainties in 
the effective area calibration. The first one was estimated by changing artificially the normalization 
of the Galactic diffuse model by ± 6 % as done in Abdo et al. (2010g). The second one is estimated 
by using modified IRFs whose effective areas bracket the nominal ones. These bracketing IRFs 
are defined by envelopes above and below the nominal energy dependence of the effective area by 
linearly connecting differences of (10%, 5%, 20%) at log 10 (£ 1 /MeV) of ( 2 , 2.75, 4), respectively. 

Using the best-fit spatial and spectral model (Dl) renormalized to the all-phase interval, the 
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Table 3. Best-fit spectral parameters obtained with gtlike using different templates for 
Puppis A above 200 MeV. The source model includes the nearby sources BckgA and BckgB 
described in Table 4. Only off-pulse photons from the Vela pulsar are selected and fluxes are not 
renormalized to the whole phase interval. The first and second errors denote statistical and 

systematic errors, respectively (see § 2.3) 


Spatial Model 

Flux(>200 MeV) 
[10 -8 ph cm -2 s -1 ] 

Photon index 

TS 

X-ray 

1.64 ± 0.23 ±0.21 

2.09 ±0.07 ±0.09 

189 

Radio 

1.64 ±0.22 ±0.22 

2.06 ±0.07 ±0.08 

187 

Uniform disk (Dl) 

1.67 ±0.23 ±0.23 

2.10 ±0.07 ±0.10 

190 

Two X-ray halves: 
Eastern half-disk 

0.68 ±0.24 ±0.17 

1.96 ±0.13 ±0.08 

82 

Western half-disk 

0.94 ±0.29 ±0.16 

2.28 ±0.14 ±0.10 

54 


Table 4. Best-fit spectral parameters obtained on the nearby sources BckgA and BckgB with 
gtlike using the uniform disk template D1 for Puppis A above 200 MeV. Only off-pulse photons 
from the Vela pulsar are selected and fluxes are not renormalized to the whole phase interval. The 
first and second errors denote statistical and systematic errors, respectively (see § 2.3) 


Source Name 

R.A. (°) 

Dec. (°) 

Flux(>200 MeV) 
[10 -8 ph cm -2 s _1 ] 

Photon index 

TS 

BckgA 

125.77 

-42.17 

1.17 ±0.26 ±0.23 

2.62 ±0.13 ±0.15 

48 

BckgB 

128.14 

-43.39 

1.72 ±0.24 ±0.20 

2.60 ±0.11 ±0.12 

89 
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7 -ray luminosity of Puppis A in the 1-100 GeV band is calculated to be 2.7xl0 34 ( D/2.2 kpc ) 2 
erg s _1 . This is among the faintest reported SNRs detected by Fermi- LAT. By comparison, SNRs 
which are interacting with dense clouds, such as IC 443 and W51C have luminosities of >10 35 erg 
s -1 (Abdo et al. 2010f, 2009a). SNR W49B is thought to have an age of ~4 kyrs, comparable to 
Puppis A, yet has a luminosity of 9.3xl0 35 erg s -1 above 1 GeV and is interacting with a dense 
molecular cloud (Abdo et al. 2010d). The relatively low GeV luminosity of Puppis A may be due 
to either less gas in the vicinity of the SNR, or to an emission mechanism different than in high 
luminosity SNRs. We explore this further in § 3 by fitting simple nonthermal emission models. 


2.4. High-Frequency Radio Spectrum 

The radio morphology of Puppis A is that of a bright shell which greatly overlaps but is not 
entirely correlated with the X-ray and IR morphologies. The integrated radio spectrum between 
19 MHz and 8.4 GHz is well fit by a single power-law with slope <a=-0.52±0.03 (Castelletti et al. 
2006). The high resolution VLA observations at 1425 and 327 MHz permitted a detailed study of 
the radio morphology and spectral index of Puppis A (Castelletti et al. 2006). While it is difficult 
to quantify, spectral flattening, to an index a ^-0.4, is observed for the radio-bright regions where 
interaction with denser gas is thought to take place. In contrast, the periphery of the SNR shell 
appears to show a steeper spectral index of a ^-0.6. The resulting global spectral index is therefore 
an average of these regions. 

We have used the 7-year all-sky data of the Wilkinson Microwave Anisotropy Probe (WMAP) 
to extend the radio spectrum of Puppis A to higher frequencies. Five bands are analyzed with 
effective central frequencies (^ e ff) of 23 to 93 GHz (Jarosik et al. 2011). To obtain WMAP flux 
densities we used template fitting of the flux-corrected 1425 MHz radio image (Castelletti et al. 
2006). This template image was smoothed with the WMAP beam profiles, and then fit in each band 
with a constant times the smoothed template plus a sloping planar baseline. To avoid emission 
from nearby Vela-X PWN, the template fit is restricted to a circular region within a 2° radius. 
The 7-year skymaps, fit template and residuals are presented in Figure 4. Fluxes are given in 
Table 5 with errors estimated from the RMS of the fit residuals in each band. For each band we 
list the effective Gaussian beam full width at half maximum (FWHM). At 23 GHz (K-band) the 
beam width is comparable to the angular diameter of the SNR. However, at higher frequencies 
the remnant is resolved as a spatially extended source. There are no previous radio detections of 
Puppis A at such high frequencies. 

The full radio spectrum of Puppis A from 19 MHz to 93 GHz is shown in Figure 5. To fit the 
radio spectrum, we exclude data below 300 MHz, which may suffer from low-frequency absorption 
by thermal electrons along the line of sight. We find a best-fit 1 GHz flux density of 141±4 Jy and 
a radio spectral index of -0.56±0.01. A radio index of -0.55 is equivalent to a nonthermal particle 
distribution index of 2.1, matching the spectral index obtained with the Fermi - LAT data in § 2.3. 
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The WMAP data above 40 GHz appear to show a decreasing flux, indicative of a spectral 
break or cutoff. Both the measured flux density at 63 GHz and the 2 a upper limit at 93 GHz 
fall below the best-fit radio spectrum. WMAP has excellent inter-band calibration (Weiland et al. 
2011), so this offers intriguing evidence for spectral curvature at high radio frequencies. To assess 
the statistical significance of the putative spectral break, we assume a synchrotron cooling break 
by A<a=-0.5, and apply the F-test to compare the y 2 fit to that of a simple power-law. We find 
a 2.8a significance for a break at 40 GHz, which is not enough to claim a detection, but highly 
suggestive. Figure 5 shows that the fit with a spectral break at ^40 GHz can reproduce the 
observed data. This is in contrast to the LAT spectrum, where there is no evidence for a break or 
cutoff at high energies. We discuss the implications of a high-frequency radio break in § 4. 


3. Nonthermal Modeling 

Determining the mechanism responsible for y-ray emission is crucial in order to measure the 
underlying relativistic particle population accelerated by the SNR. The excellent correlation be- 
tween the GeV and X-ray morphologies is suggestive of potential y-ray emission mechanisms. First, 
the shock heated thermal X-ray gas produces strong IR emission by heating dust grains (Arendt 
et al. 2010). Relativistic electrons that inverse Compton (IC) scatter off this IR emission, as well 
as off the cosmic microwave background can contribute to y-ray emission. Alternatively, the in- 
creasing X-ray brightness from the South-West to the North-East of the remnant is explained by 
an increasing density of swept-up gas, rising from 0.5 to 4 cm -3 . The correlation with gas density 
suggests either nonthermal Bremsstrahlung emission by relativistic electrons or neutral pion decay 
by relativistic protons, as both scale with increasing density. 

Here we adopt the simplest possible assumption that all emission originates from a single 
population of accelerated protons and electrons contained in a region characterized by a constant 
matter density and magnetic field strength. The emitting region is assumed to be homogeneous 
with a volume given by V — /(47t/3) i? 3 where / < 1 denotes a filling factor and R — 15 pc is 
the radius of the remnant. For simplicity, the filling factor is fixed at / = 1 in the following. The 
particle spectra are assumed to follow a power-law with an exponential cutoff dN/dE oc r\ e ^ E _r 
x exp(-E/E max ), with the same spectral index and energy cutoff for both electrons and protons. 
Electrons suffer energy losses due to ionization, Bremsstrahlung, synchrotron processes and IC 
scattering. The modification of the electron spectral distribution due to such losses was calculated 
according to Atoyan et al. (1995), where electrons are assumed to be injected at t = 0 from an 
impulsive source. The particle spectrum is then evolved for 3700 years, the age of Puppis A. 

We consider three models in which each of the three plausible emission mechanisms (IC scatter- 
ing, Bremsstrahlung, proton-proton interaction) dominates. In fitting the LAT spectrum of Puppis 
A, the normalization, maximum energy cutoff and magnetic field are left as free parameters. Since 
we assume that protons and electrons have identical injection spectra, the spectral index of the 
particle spectrum below the cutoff energy is determined and fixed at ~ 2.1 by modeling the radio 
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spectrum as synchrotron radiation by relativistic electrons. The spectral break in the WMAP data 
above 40 GHz being marginal with the current statistics, we do not try to reproduce this break 
in the following scenarios proposed. Modeled spectral energy distributions (SEDs) are presented 
in Figure 6. Model parameters are given in Table 6. We also give the total energy of accelerated 
particles integrated above 1 GeV for protons, and above 511 keV for electrons. While the chosen 
parameters are not unique in their ability to fit the broadband spectrum, they are representative 
of the difficulties and advantages of each scenario. We discuss the viability of each emission model 
in the following subsections. 


3.1. Inverse Compton Dominated Model 

Puppis A has several photon fields which can be up-scattered by relativistic electrons to pro- 
duce IC emission: the cosmic microwave background radiation (CMBR), the infrared continuum 
produced by the collisional excitation of swept-up dust, and the local interstellar radiation field 
(ISRF). The luminosity of IC y-rays is proportional to n p h e 1 / 2 , where n^ is the average photon 
density and e is the average photon energy (Gaisser et al. 1998). Typically, the CMBR is the 
dominant radiation field, characterized by n p h ~ 400 cm -3 and e ~ 6.3 xlO -4 eV. Dust emission 
from Puppis A is fit by a two temperature black body with Ti = 150 K, n p h—l2 cm -3 , e=0.04 eV; 
and T2=45 K, n p h—20 cm -3 , e=0.01 eV (Arendt et al. 1991). Summing these components, the 
total contribution to IC y-rays from IR dust emission is only ^60% that from the CMBR. As Pup- 
pis A lies at a Galactocentric radius comparable to the Sun, we assume the local ISRF spectrum 
(Mathis et al. 1983) which peaks at 1.2 eV, and contributes ^25% compared to the CMBR. While 
shocked dust and the stellar optical background contribute significantly to IC emission, they do 
not dominate over the CMBR, and cannot solely explain emission from Puppis A. 

Accounting for these photon fields, an IC-dominated model with a particle index of 2.1 and a 
cutoff at 500 GeV provides a good fit to the data. With a magnetic field of ^8 fi G, this scenario 
requires a total energy in CR electrons of 2.9xl0 49 erg. It should be noted that the electron-to- 
proton ratio is not well constrained and was fixed here at r] e /rj p = 1. However, this value needs 
to be larger than 0.1, and thus in excess of the ratio found for local cosmic-ray abundances, to 
inject a reasonable energy content in radiating electrons. In the same way, the average gas density 
needs to be lower than 0.3 cm -3 to reduce the Bremsstrahlung component. This limit is very close 
to the lower limit from thermal X-ray observations (Hwang et al. 2005) and is somewhat low for 
unshocked ISM in the vicinity of a molecular cloud. 


3.2. Bremsstrahlung Dominated Model 

Bremsstrahlung emission, from CR electrons encountering gas of moderate density, can provide 
a reasonable emission mechanism. X-ray and IR observations constrain the density in the shock 


-13- 


Table 5. WMAP Flux Density for Puppis A. Five bands are analyzed with effective central 

frequencies (z/ e ff) of 23 to 93 GHz (see § 2.4). 


Band 

v eff 

(GHz) 

Flux Density 

(Jy) 

FWHM 

(°) 

K 

22.7 

25.8±2.3 

0.93 

Ka 

33.0 

20.8±1.7 

0.68 

Q 

40.6 

17.8±1.5 

0.53 

Y 

60.5 

12.3±1.5 

0.35 

W 

93.0 

3.0±2.3 a 

0.23 


a The source is not detected at 93 GHz. We 
find a 2 a upper limit of 7.6 Jy. 


Table 6. One-Zone Model Parameters 


Model 

Index 

-^max 

[TeV] 

n H 

[cm -3 ] 

B tot 

[mG] 

Ve/Vp 

W p 

[erg] 

W e 

[erg] 

IC 

2.1 

0.5 

0.3 

8 

1 

8.0xl0 48 

2.9xl0 49 

Brems. 

2.1 

0.5 

4 

13 

1 

3.5 xlO 48 

1.3xl0 49 

7r°-decay 

2.1 

0.8 

4 

35 

0.02 

4.0xl0 49 

2.8xl0 48 
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interaction region to be ^4 cm -3 (Arendt et al. 2010). Assuming an electron-to-proton ratio of 
1 and a gas density of 4 cm -3 , the best-fit Bremsstrahlung dominated model has a particle index 
of 2.1 and a cutoff above 500 GeV. The total energy in CR electrons is 1.3xl0 49 erg, and the 
magnetic field is ^13 fi G. Given that the SNR is adiabatically expanding, compressing the swept- 
up gas and magnetic field lines by a factor of ^4, a magnetic field of 13 /iG is reasonable for 
compression alone in the post-shock region from which radio and y-ray emission may arise. We 
also note that at densities of ^4 cm -3 , Bremsstrahlung emission dominates over hadronic 7r° decay 
for electron-to-proton ratios greater than ^0.1, again in excess of cosmic-ray abundances. 


3.3. Hadron Dominated Model 

We then consider a 7r°-decay model to account for the broadband spectrum of Puppis A. 
The flux and spectrum of the y-rays produced by 7r°-decay are calculated following the analytical 
approximations by Kelner et al. (2006). Here, we assume a gas density of 4 cm -3 , and an electron- 
to-proton ratio of 0.02, consistent with the locally observed cosmic-ray composition. The particle 
index is fixed at 2.1 with a cutoff above 800 GeV. The total energy in CR hadrons is 4.0xl0 49 erg, 
and for electrons is 2.8 xlO 48 erg. Assuming a typical explosion energy for Puppis A, this hadronic 
model indicates that a few percent of the initial kinetic energy is transferred to CRs within a few 
thousand years. The magnetic field value required to reproduce the synchrotron radio observations 
is ^35 /iG, assuming that the radio and y-ray emissions arise from the same volume. However, if 
the y-ray emission originates from higher density gas with a lower volume filling fraction than the 
radio emission, the average magnetic field would be less than our fitted value. 

From the three emission models above, it is clear that leptonic scenarios require an electron-to- 
proton ratio at least ten times larger than the standard cosmic-ray abundance ratio. In addition, 
the inverse-Compton dominated model requires an extremely low density as well as a relatively low 
magnetic field value. However, with the current statistics, it is not possible to discriminate between 
these different models, and the parameters listed in Table 6 provide an estimate of the energetic 
particle population in the SNR Puppis A under reasonable assumptions. It should be noted that 
the CR content (electrons and protons) estimated in the leptonic models should be taken as a lower 
limit since the electron-to-proton ratio is fixed at a maximum value of y e /yp = 1 in the fit. 


4. Discussion 

Puppis A is adiabatically expanding into the interstellar medium in the vicinity of a large 
molecular cloud. At an age of ^3700 years, the SNR has not yet encountered the bulk of the large 
cloud, but is likely to at some point in the future. This makes Puppis A an interesting transitional 
case between young SNRs still evolving into a circumstellar medium (e.g. Cas A), and older SNRs 
which are interacting with large, dense molecular clouds (e.g. IC 443). At an age of only 350 years, 
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Cas A is detected at both GeV and TeV energies with a 7-ray spectral index of ^2.1 (Abdo et al. 
2010c; Acciari et al. 2010). Multi-zone modeling of the SED found that a hadronic CR content 
of at least 6xl0 49 erg is required (Araya & Cui 2010). This is comparable to the total energy in 
nonthermal protons in Puppis A estimated in our hadronic model. In contrast, SNR IC 443 has a 
luminosity one order of magnitude greater than Puppis A, with a total energy content in protons 
of only ^few times 10 49 erg (in hadronic models), similar to Puppis A. The difference in luminosity 
can be easily explained by the much higher average density in IC 443, as it is interacting with large 
molecular clouds. Cas A is thought to have a target density comparable to that for Puppis A, and 
so a comparable luminosity is to be expected. 

It is interesting to speculate on how future observations may constrain the nature of the 7-ray 
emission in Puppis A: either with improved statistics from Fermi - LAT or from a detection at TeV 
energies. Extrapolating from the LAT flux and index, the uniform disk has an estimated photon 
flux of 2 xlO 11 ph cm 2 s 1 above 0.1 TeV, equivalent to ^5% of the Crab flux. The current 
generation of Cherenkov telescopes should be capable of detecting Puppis A, if no high energy 
cutoff is present. If TeV emission is detected toward Puppis A, we expect it to be brightest in the 
East where there is also enhanced radio and X-ray emission from interaction with higher density 
ISM toward the adjacent molecular cloud. As can be seen in Figure 2, the centroid of LAT emission 
moves toward the East with increasing energy. The differences in 7-ray photon index agree with 
variations in the observed radio spectral index (Castelletti et al. 2006). The bright Eastern shell 
appears flatter (a 0.5, T ^2), while the fainter Western extent is generally steeper (a 0.T, 
T ^2.4). TeV observations can provide an important test 

Future radio observations are also needed to confirm a spectral break or cutoff at >40 GHz 
indicated by the WMAP spectrum. Few SNRs have been convincingly shown to have breaks at 
high frequencies. One such SNR is S147, which has a break at ^1.5 GHz. It was also recently 
identified as a spatially extended Fermi - LAT source (Katsuta et al. 2012). To explain the radio 
and 7 ray spectra, a two-zone model was invoked: diffuse, low density gas dominates the radio 
emission at low frequencies, while high-density shock-compressed filaments give rise to 7-rays and 
radio emission at frequencies above the radio break. 

However, there are significant differences between S147 and Puppis A. S147 is an order of 
magnitude older and has entered the radiative phase (indicated by dense H a filaments), while 
Puppis A is still non-radiative, except for a few bright X-ray knots. It is also more difficult to explain 
the break frequency in Puppis A. For a magnetic field of an electron radiating synchrotron 

emission at a peak frequency z/ghz has an energy E = 14.7 (^ghz/B^g ) 0 ' 5 GeV (Reynolds 2008). 
The radio break in S147 occurs at ~1.5 GHz, and may be explained by synchrotron cooling over 
the long lifetime of the remnant. It could also be the effect of damping of magnetohydrodynamic 
turbulence (due to ion-neutral collisions) in a model where shock re-accelerated CR particles are 
adiabatically compressed and energized as a result of radiative cooling behind the cloud shock 
(Maikov et al. 2011). No radio break is observed for the denser optical filaments. A synchrotron 
break at ^40 GHz would require a magnetic field in excess of 1 mG to have existed for the entire 
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lifetime of Puppis A, in order to cool the electron spectrum. However, no corresponding break is 
observed in the LAT spectrum with the current statistics, so two electron populations would be 
required to fit the radio and 7-ray data. 

Clearly, further observations to characterize the radio and 7-ray spectra are needed. High- 
frequency radio observations are needed to spatially resolve the radio spectral break or cutoff. 
Increased sensitivity with the continued observations by Fermi - LAT will allow any differences 
between the Eastern and Western regions to be firmly established. Observations at TeV energies 
of Puppis A will be interesting in several respects: (1) to determine whether a spectral cutoff or 
break is present at higher energies and perhaps discriminate between emission mechanisms, (2) to 
measure differences between East and West of the SNR, and (3) to search for escaping CRs which 
may have encountered the nearby Eastern cloud. 


5. Conclusions 

We report on high energy 7-ray emission from the Puppis A SNR detected with Fermi- LAT. 
The source is clearly spatially extended at energies > 800 MeV. The morphology of GeV emission is 
well correlated with X-ray and IR morphologies. Spectrally, the SNR is well described by a power- 
law index of 2.1 in the LAT energy range. We also report the detection of Puppis A with WMAP at 
20 to 93 GHz. Extending the radio spectrum to high energies reveals a putative spectral break or 
cutoff at a frequency above 40 GHz, which would have interesting implications for the nonthermal 
electron population if confirmed. We fit the radio to 7-ray SEDs using IC, Bremsstrahlung and 
hadronic dominated models. All emission mechanisms are able to fit the data, though with different 
magnetic field strengths and energetics of relativistic particles. For the hadronic models a total 
energy in CR protons of ^4xl0 49 erg is needed, whereas for leptonic models, at least (1 - 3)xl0 49 
erg in relativistic electrons is required. While Puppis A has a relatively low 7-ray luminosity in 
comparison to other SNRs identified by Fermi , this may be explained by the low average gas 
density, nn ^4 cm -3 , swept-up by the SNR. Future observations at high radio frequencies and at 
higher 7-ray energies will help to differentiate between leptonic and hadronic emission models. 
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2 7 15 27 42 *0 83 107 13G 

Fig. 1. — Test Statistic (TS) map in celestial coordinates of the SNR Puppis A. Only photons 
with energies above 800 MeV in the off-pulse window of the Vela pulsar are selected. The TS 
was evaluated by placing a point-source at the center of each pixel, Galactic diffuse emission and 
nearby sources being included in the background model, except the two background sources BckgA 
and BckgB which are indicated with yellow crosses. The location of the compact object inside 
Puppis A, PSR J0821— 4300, is indicated by a black diamond. The gray circle corresponds to the 
best fit model (Dl) provided by point like in Table 1, the dashed line showing its division in two 
half-disks along the North/South line in celestial coordinates. The color-coding is represented on 
a square-root scale. 
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Fig. 2. — Test Statistic (TS) map of the SNR Puppis A above 800 MeV in the off-pulse window of 
the Vela pulsar (Top) and above 5 GeV in the all-phase interval (Bottom). The TS was evaluated 
by placing a point-source at the center of each pixel, Galactic diffuse emission and nearby sources 
being included in the background model. The location of the background source BckgA added to 
the source model is indicated with a yellow cross. Green and cyan contours correspond to images 
at different wavelengths: ROSAT HRI X-ray emission contours (Green, Left) and 1.4 GHz radio 
continuum contours (Cyan, Middle). The contour levels are 0.02, 0.19, 0.37, 0.55 counts s -1 m 2 
for X-ray emission, and 5, 10, 15, 20 mJy beam -1 for 1.4 GHz radio continuum. The yellow disk 
corresponds to the best fit model (Dl) provided by point like in Table 1. The color-coding is 
represented on a linear scale for all images. 
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Fig. 3. — Fermi - LAT spectrum of Puppis A using the X-ray spatial model described in the text. 
The bold-dashed (red) line shows the fit of a power law to the overall spectrum derived from all of 
the data with energy E > 200 MeV. The data points (crosses) indicate the flux measured in each 
of the seven energy bins renormalized to the whole phase interval, as indicated by the extent of 
their horizontal lines. The statistical la errors are shown in dark (blue), while the light (red) lines 
take into account both the statistical and systematic errors as discussed in § 2.3. The gray shaded 
area shows the extent of the Fermi 68% confidence band for the whole SNR; the two dashed (blue) 
areas show the extent of the Fermi 68% confidence band for each half-disk of the SNR. 
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Fig. 4. — Template fits to the WMAP data within a 2° radius of Puppis A. Five columns are labeled 
with the corresponding WMAP radio bands. Three rows present the 7 year skymap images, models 
resulting from a fit of the 1.4 GHz radio template plus a sloping planar baseline, and the relative 
residuals (defined as the fit residual divided by the model). The upper and lower limits of the linear 
colorbar are given beneath each image. Data and models are on the same scale. Puppis A is clearly 
detected in four bands, but not in the highest frequency W-band (93 GHz). We note that the radio 
point source to the north of Puppis A, which is clearly visible in the residuals, corresponds to the 
massive star forming region, IRAS 08212-4159. Including this source in our template fit does not 
affect the fitted values for the SNR. All images are in celestial coordinates. 
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Fig. 5. — Radio spectrum of Puppis A from the flux densities listed in Table 5 and values compiled 
by Castelletti et al. (2006). The solid line shows the spectrum for the entire SNR with a flux density 
of 130 Jy at 1 GHz and a spectral index of <a=-0.56. The filled triangle shows the 2a upper limit 
at 93 GHz. The inclusion of WMAP data shows a clear spectral break for frequencies above ^40 
GHz. The dashed line shows a spectral break at 40 GHz from an index of -0.55 to -1.05 which is 
consistent with the WMAP data. 
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Fig. 6. — SED models of both the Fermi - LAT spectrum in Figure 3 and the radio spectrum in 
Figure 5, for which IC (top), Bremsstrahlung (middle) and 7r°-decay (bottom) are the dominant 
emission mechanism (see Table 6). In each model the radio data are fit with a synchrotron compo- 
nent. All models show the contributions of 7r°-decay (long dashed, red), Bremsstrahlung (dotted, 
blue), and IC emission (dashed, blue) from CMB, IR dust photon field, and stellar optical photons. 
The sum of the three y-ray components is shown as a solid black curve. All data points are shown 
with 1 a error bars. The downward arrow shows the 2 a upper limit from WMAP at 93 GHz. 


